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If the hghter top squark has no kinematically ahowed two-body decays that con- 
serve flavor, then it wifl form hadronic bound states. This is required in models that 
are motivated by the supersymmetric Uttle hierarchy problem and obtain the correct 
thermal relic abundance of dark matter by top-squark-mediated neutralino annihi- 
lations, or by top-squark- neutralino co-annihilations. It is also found in models that 
can accommodate electroweak-scale baryogenesis within minimal super symmetry. I 
study the prospects for detecting scalar stoponium from its diphoton decay mode at 
the Large Hadron Collider, updating and correcting previous work. Under favorable 
circumstances, this signal will be observable over background, enabling a uniquely 
precise measurement of the superpartner masses through a narrow peak in the dipho- 
ton invariant mass spectrum, limited by statistics and electromagnetic calorimeter 
resolutions. 
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I. INTRODUCTION 

The classic collider signatures for supersymmetry depend on the presence of missing energy 
carried away by a stable, neutral, weakly interacting lightest super symmetric particle in each 
event. (For reviews, see [H, 0, S]- This paper follows the conventions and notations of the last 
reference.) In a sense, this is disappointing, because there are no kinematic mass peaks whose 
positions would yield measurements of superpartner masses. Observables at hadron colliders can 
give precision determinations of superpartner mass differences by way of kinematic edges and other 
distributions, but if i?-parity is conserved the overall mass scale will be much harder to ascertain 
with precision in most models 0, 0] • 

A possible exception occurs if the super symmetric particles can form resonances that annihilate 
into final states containing only Standard Model particles with strong or electromagnetic interac- 
tions. An example is stoponium, a bound state of a pair of top squarks. The lighter top squark (or 
stop) tl has possible flavor-preserving two-body decays into the lightest chargino or neutralino, 

h ^ bCu (1.1) 

ii ^ m, (1.2) 

which, if kinematically allowed, would cause it to decay long before it could form a hadronic bound 
state. However, the first decay will be closed if the chargino Ci is not at least 5 GeV lighter 
than the top squark, and the large top quark mass means that the second decay may well also 
be kinematically closed. In most of the so-called mSUGRA parameter space, this situation is not 
encountered, but it is nevertheless quite possible and even common in other model frameworks. 
Then one must consider three-body (and four-body) decays that preserve flavor, and a two-body 
decay to charm that violates flavor: 

ii W''*hNi, (1.3) 
ii cNi. (1.4) 

The partial widths associated with the decays (II. 3p and (jl.4p are known 0, [tI to be far smaller 
than the binding energy of stoponium or other bound hadronic states involving ti. 

Therefore, it is worthwhile on general grounds to consider signals for the production and decay 
of stoponium. At hadron colliders, stoponium is produced primarily in gluon-gluon fusion, with 
the largest cross-section for the IS* {J^'~' = 0~^'^) scalar ground state, denoted in the following as rj^. 
This state will decay primarily by annihilation, with the possible two-body final states including 
gg, 77, W~^W~, ZZ, h^hP, ti, bb, and NiNi. The QCD backgrounds for the gluon and quark 
final states are too huge to contemplate a signal. Also, the W~^W~ , ZZ, and hPhP final states are 
plagued by either large backgrounds or small branching fractions, and the NiNi final state does 
not give a reconstructable signature. However, Drees and Nojiri in [1, 0] pointed out that rfi 77 
provides a viable signal at the CERN Large Hadron Collider. (See also refs. flo|-[l3| for earlier 
works related to stoponium at hadron colliders.) 

Refs. [1, 0] appeared before two important experimental developments which bring this possi- 
bility into sharper focus. The first is the 1995 discovery of the top quark. Second, the results of 
WMAP and other experiments have bracketed the density of cold dark matter in the standard cos- 
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mology This is important because the measurement of JIdm^^ ~ 0.11 can be correlated 

with a top squark light enough to forbid the decays (jl.ip and (|1.2p in at least two scenarios. First, 
neutralino LSP and top-squark co-annihilations c an g ive the observed dark matter density if the 



mass difference m^^ — m^y^ is in a small range 18|- 2^. Second, the neutralino LSP can efficiently 



pair-annihilate into a top-anti-top pair, mediated by the i-channel exchange of a top squark that 



is not more than about 100 GeV heavier than the LSP [231 . l24f |. These two scenarios are often 
continuously connected in parameter space, but the latter one requires far less fine adjustment of 
parameters to realize. It does, however, require that the gaugino masses are not unified in the 
way assumed in mSUGRA models. A particularly attractive model framework that meets the re- 
quirements of the stop-mediated annihilation scenario is provided by "compressed supersymmetry" 



231 . l24f | , in which the gluino mass parameter is taken to be significantly smaller than the wino mass 
parameter at the scale of apparent gauge coupling unification Mgut ~ 2 x 10^^ GeV. A reduction 
of the gluino mass compared to the wino and bino mass parameters can also ameliorate |25] the 
supersymmetric little hierarchy problem. Another quite different motivation for a light top squark 
is provided by models that can have a strongly first-order phase transition leading to electroweak- 



scale baryogenesis 26|- 28|]; this can also incorporate the neutralino-stop co-annihilation scenario 
for dark matter. 

In this paper, I will consider the 7]^ 77 signal at the LHC along similar lines to ref. 0], taking 
into account the now known top-quark mass, considering motivated models that agree with the 
observed dark matter density and Higgs mass constraints, and using a more liberal angular cut 
but a more conservative energy resolution for the electromagnetic calorimeter. I also correct (see 
Appendix) factor of 2 errors appearing in the gg and 77 partial decay widths in ref. 0|. This 
leads to a somewhat more pessimistic evaluation of the detection potential, but in many parts of 
parameter space the diphoton signal will be detectable given a large integrated luminosity at the 
LHC. 



II. DIPHOTON SIGNAL FOR STOPONIUM AND BACKGROUNDS 

The leading-order partial decay widths of rj^ into gluon and photon final states are: 

nvi^gg) = ^a||i?(o)|V<-, (2.1) 

nVi^ll) = ^«2|i?(0)|Vm2., (2.2) 
where i?(0) = ^/AiripiO) is the radial wavefunction at the origin. In the Coulomb app roximation 



to the bound state problem, |i?(0)p/m^_ = 4a|m^-/27. However, the study of ref. [29[ indicates a 



softer potential, with the Coulomb limit not obtained even for very large bound state masses. In 
the following, I will adopt the = 300 MeV parameterizations of the wavefunction at the origin 



and the binding energy as given in ref. 29l |: 



|i?(0)|V"i^j = (0.1290 + 0.0754L + 0.0199L2 + 0.0010L3) GeV, (2.3) 
2m^^ - m^- = (3.274 + 1.777L + 0.560^^ + O.OSIL^) GeV, (2.4) 



4 



FIG. 1: The binding energies for the IS, 2S, 
and IP stoponium states as a function of 
the stoponium mass, as computed from the 
potential model of ref. with = 300 
MeV. 
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where L = ln(m^^/250 GeV). The binding energies of the IS ground state and the 2S and IP ex- 
cited states are shown in Figure [H It should be noted that these results are based on a considerable 
extrapolation from known experimental results on cc and bb bound states, and other potentials can 
give quite different results. (For example, see the ones reviewed in ref. 0].) The partial width 
into gluons is of order 2 MeV over the considered range of m^,. The fact that the binding energy 
is much larger shows that the stoponium bound state will indeed form, provided that other partial 
widths do not overwhelm r(7y^ — > gg) by a factor of 1000, a requirement easily satisfied by models 
studied below. 

At leading order and in the narrow-width approximation, the production cross-section for rjf; in 
pp collisions is given in terms of its gluonic decay width by 



a{pp rii) = -^^^im 99) J dx-g{x, Q^)g{T/x, Q^), (2.5) 



where g{x, Q^) is the gluon parton distribution function, and r = tci^J s in terms of the pp collision 



energy squared s. In the following, I use the CTEQ5L 3l| set for the parton distribution functions, 
evaluated at Q = ?ti^- for the signal and Q = M^^ for the backgrounds. 

For comparison, the ratio of the stoponium production cross-section to that of a Standard Model 
Higgs boson H with the same mass is just a{pp — > rji)/a{pp H) = T{r]^ gg)/V{H — > gg) at 
leading order. For m^- = (200, 250, 300, 400, 500, 600, 700, 800) GeV, this is approximately a{pp 
Vi)/cr{pp ^ H) = (1.44,0.74,0.40,0.098,0.064,0.054,0.050,0.048). However, BR(r/£ 77) is much 
larger than BR(f/' 77) for masses larger than 140 GeV, because the latter is loop-suppressed 
compared to relatively huge widths into ZZ^*\ WW^*\ and tt final states. This explains why the 
stoponium signal r/^ 77 can be viable over the mass range where if ^ 77 observation is not 
possible. In contrast, BR(?7|- — > ZZ)/BK{H — > ZZ) turns out to be at most about 0.3 over the 
same mass range. This explains why the search for stoponium in — > r/^ — > ZZ in the 
and i^i^uu channels is almost certainly not viable for masses where it is the best search option 
for H. 

Now, multiplying eq. ()2.5p by the branching fraction into a diphoton state, and rearranging the 
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factors, one obtains: 



a{pp r]i 77) = ^— 3-BR(t/^- gg)T{r]i 77) j dx-g{x, Q'^)g{T/x, Q^). (2.6) 



This way of writing the result is useful because, in many realistic models, the gluonic decay dom- 
inates over all other final states. Therefore, it is instructive to adopt an idealized limit where 
BR(r/j gg) ~ 1 as a standard reference scenario. (Note that the 77 partial width is typically 
about 0.005 of the gg partial width.) Then results for particular models can be obtained by scaling 
the signal cross-section by the actual BR(r/^ — > gg). 

We next consider the diphoton backgrounds at the LHC. The pp 77 process has parton-level 
contributions: 

qq 77, (2.7) 
gg 77, (2-8) 



with leading-order differential cross-sections found in [3J]- 35|]. The leading order total cross-section 
for eq. (j2.7p is proportional to ln[(l -|- zo)/{l — zq)] — zq, where zq is the cut on | cos 6*^,1, with 6^ the 
photon momentum angle with respect to the beam direction in the center-of-momentum frame. 
Since the signal is isotropic, with a total signal cross-section proportional to zq, it follows that 
S/\fB is maximized for zq « 0.705. The process (12. Sh involves Feynman diagrams that have quark 
box loops. It is somewhat more central than the qq background, and so would favor a larger cut 
Zq, but it is quite subdominant over the diphoton mass range considered here. For simplicity, I will 
impose a cut on both signal and background of 

I cos 6** I < 0.7 (2.9) 

in the center-of-momentum frame. This guarantees a high px for the photons, for high-mass 
stoponium states. In addition, I will require that the photons be well-separated from the beam 
direction and the remnant beam jets, so 

|cos6'|<0.95 (2.10) 

(or |r/| < 1.83) in the lab frame. The cut (j2.9p is apposite for low rapidities, and (j2.10p for high 
rapidities. The results for the background at leading order after these cuts are shown in fig. [2] as 
a function of the invariant mass of the diphoton system, M^^. 

Higher-order corrections to the background can be quite important, increasing the cross-section 



after cuts by a factor of 2 or more |36[|-[39[|. This includes a large contribution from the hard scatter- 
ing process qg — > 775 (and qg — > 77^) and the related process qg — > ^q followed by a photon from 
the fragmentation of the quark jet, as well as from double fragmentation contributions. Imposing 
isolation cuts on hadronic activity near the photons and requiring the absence of additional hard 
jets reduces these backgrounds considerably. The isolation cut requirement is necessary anyway, 
to eliminate an otherwise large background from jets faking photons, including vr'' and ry decays 
that are not resolved in the electromagnetic calorimeter. Also, the initial-state gluon contribu- 
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FIG. 2: The differential cross-section 
da/dM^j for the diphoton backgrounds in 
pp collisions at y/s = 14 TeV, due to 
the parton-level processes qq 77 and 
99 ~* 77 1 8.t leading order. Here Mj^ is the 
diphoton invariant mass. The cuts imposed 
on the angle with respect to the beam axis 
are | cos0*| < 0.7 in the diphoton center-of- 
momentum frame and | cos9\ < 0.95 in the 
lab frame. The kink in the gg background 
at M-y-y — 2mt is due to the threshold in 
the top-quark box loop. 
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tions are not as important in the present case as in the well-studied Standard Model Higgs signal 
mass range mn < 140 GeV, because the gluon parton distribution function is relatively smaller 
at larger x. The background contributions from jets without charged tracks faking photons can 
likely be reduced to a small level with isolation cuts imposed in offline analysis 40]. In any case, 
because the stoponium diphoton resonance will be very narrow, in practice the background should 
be determined directly from LHC data by a sideband analysis. 

Inclusion of the higher-order background contributions is beyond the scope of this paper, es- 
pecially since the corresponding higher-order corrections to the stoponium signal cross-section are 
not available. In addition, the increased background will likely be at least partly compensated 
for by contributions to the signal from production of excited stoponium states, followed either by 
decays to the IS r/^ state by emission of photons or soft mesons or by direct decays to 77 For 
example, the 25 state has a binding energy that is probably only slightly less than the IP state (see 
fig- H)) so decays of 2S stoponium to the IP state and a meson would be kinematically forbidden 
according to this potential model. The 2S non-annihilation decays therefore will most likely go 
entirely to the IS" ground state. These signal contributions will be effectively merged due to the 
detector energy resolution, leading to an overall enhancement of the signal of perhaps a factor of 
1.5 (see fig. 9 of ref. [^). This unknown enhancement is not included here, to be conservative. 

The energy resolutions of the ATLAS and CMS electromagnetic calorimeters at the LHC will 
clearly dominate over the very small intrinsic width of stoponium in determining the experimental 
width of the signal peak. Therefore, to estimate the significance of the signal over the background, I 
consider a bin with width 0.04m^-, chosen to contain essentially all of the signal peak on which it is 
centered. (See ref. [40|] for a CMS study of a similar diphoton signal peak at higher masses, and the 
CMS and ATLAS physics technical design reports 0, 0] for estimates that put the electromagnetic 
calorimeter resolutions at roughly the per cent level.) The resulting comparison of the signal and 
background within such a bin is found in fig. [3l Here I have taken the signal corresponding to the 
idealized reference case of BR(r/j — > gg) + BR(r/j 77) = 1. 

The resulting integrated luminosities needed to reach significances S/\fB = 2,3,4,5, taking 
enough events to use Gaussian statistics, are shown in fig. [H Because there are many bin widths 
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FIG. 3: The pp rj^ 77 cross-section 
and the irreducible background in a bin 
with \Mjj — m,,j| < 0.02m^-, for pp col- 1q' 
lisions at ^/s = 14 TeV, as a function of ' ' 
the stoponium mass m^-. Both are com- 
puted at leading order, with the same angu- 
lar cuts as in Figure [21 The signal assumes 
an idealized limit in which BR(77j gg) + 
BR{r]t 77) is 100%. 
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FIG. 4: Total integrated luminosity yield- 
ing an expected S/\/B — 2,3,4,5 for 
in a bin with \M^^ — m,,j| < 0.02to,,j, as a 
function of the stoponium mass m^^j, for pp 
collisions at ^/s = 14 TeV, in the idealized 
limit that BR{r]f — > gg) + BR(7yj 77) is 
100%. The integrated luminosity needed to 
achieve a given S/^/B can be obtained by 
scaling by l/[BR{r]t gg)]^. 
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over the mass range considered, the probability of a 2-sigma excess in one of them arising just 
from a fluctuation in the background is not negligible. However, the LHC will likely have already 
produced a preliminary estimate of the ii mass from gluino or direct open stop production before 
the stoponium diphoton signal becomes feasible, so the search range for the stoponium mass peak 
will not be too large. The luminosities in fig. H] should be multiplied by l/[BR(r/£ 99)]'^^ since 
the required luminosity scales like the square of the signal cross-section eq. (12. 6p . 

From fig. [4] one sees that the expected significance for a 500 GeV stoponium resonance will be 
at most only about S/\^ = 2 for a canonical high luminosity year of data (100 fb~^). This is 
more pessimistic than in ref. |9|], due in part to a factor of 2 error in that paper in the rj^ — > 77 
width, but also due to their assumption of electromagnetic calorimeter resolution providing an 
acceptable bin for M^^ that is twice as narrow as assumed in the present paper. However, a more 
sophisticated approach based on a maximum likelihood fit, which is beyond the scope of this paper, 
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will certainly do better than the simple counting in a single bin used here. (Note that the detector 
mass resolution will in any case be smaller than the bin width needed to catch all of the signal 
events.) Also, the integrated luminosity needed is proportional to the square of the signal cross- 
section and to the reciprocal of the background, and so is strongly dependent on assumptions that 
are difficult to evaluate confidently at present. 



III. RESULTS FOR COMPRESSED SUPERSYMMETRY MODELS 

In the previous section, I estimated the integrated luminosity needed to achieve detection of the 
stoponium resonance at a given significance, but considering an idealized reference model where 
BR(r/^ gg) was nearly 100%. In this section and the next, I consider the actual branching 
ratio achieved in realistic, motivated models that satisfy the dark matter density and Higgs mass 
constraints. As discussed in the previous section, the integrated luminosity needed for discovery 
scales like l/[BR(?7j- gg)]'^ 



First, I will follow ref. 23|, l2j] and consider models where the bino, wino, and gluino masses 



can be parameterized at Mgut by: 

Ml = mi/2(l + C24), (3.1) 

M2 = mi/2(l + 3C24), (3.2) 

Ms = mi/2(l - 2C24), (3.3) 

corresponding to an F-term source for supersymmetry breaking in a linear combination of the 
singlet and adjoint representations of SU{f>). For the sake of simplicity, I also assume a common 
scalar mass mo and scalar trilinear coupling ^0 ^-t Mgut- For C24 of order 0.2, one finds that 
the supersymmetric little hierarchy problem is ameliorated, with a significant part of parameter 
space where ti is the next-to-lightest supersymmetric particle. The dark matter thermal relic 
abundance can be sufficiently suppressed by NiNi — > tt due to t-channel ti exchange, giving 
^DM^^ = 0.11 in accord with observations. I impose this as a requirement, by adjusting the 



value of rriQ for fixed values of the other parameters, using the program micrOMEGAs 2.0.1 41 1 



(checked for approximate agreement with DarkSUSY [42]) interfaced to the supersymmetry model 



parameters program SOFTSUSY 2.0.11 43| (checked for approximate agreement with SuSpect 44 1 
and IS A JET {451). The resulting mo values are reasonably small and do not require fine-tuning. In 
these models, m^^ — m^^ < 100 GeV, so that stoponium will indeed form as a bound state before 
it has a chance to decay. ^ These stop- mediated annihilation models are continuously connected 
in parameter space to more fine-tuned regions in which the ti, A'^i mass difference is just right to 
allow efficient stop-neutralino co- annihilations. 

In general, as shown in ref. fo^, the most important final states in competition with the gg 
and 77 ones are W^W~ , ZZ, h^h^, tt, bb, and NiNi. Formulas for the decay widths for these 
final states were given in ref. 0], and are presented in the Appendix of the present paper in a 



t The LHC phenomenology (other than from stoponium) and dark matter detection prospects of these models have 
been discussed in refs. [j^. [23. [i^ . There has recently been considerable interest in the phenomenology of other 
models that achieve realistic dark matter with non-universal gaugino masses [47t-[63l]. 
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FIG. 5: The branching ratios of scalar stoponium into the most important final states, for some representative 
models of the type described in the text. In all cases, tan/? — 10, /i > 0, and M\ varies, with mo adjusted 
to give f2DM^^ — 0.11. The upper left, upper right, lower left, and lower right panels have respectively 
(C24, -Ao/Mi) (0.19, 1), (0.21, 1), (0.21, 1.5), and (0.24, 1.5). The thicker part of the gg line indicates the 
range of stoponium mass for which stop-mediated annihilations N\N\ —* tt contribute more than 50% to 

1/^^DM^^- 



different notation. Results for the branching ratios in four typical model lines are shown in fig. 
[5] These model lines each have a continuously varying overall gaugino mass scale, parameterized 
by the bino mass parameter Mi, with the wino and gluino masses at Mgut then given by fixed 
values of C24, as in eqs. (I3.ip - (l3.3p . The parameters tan/3 = 10 and Aq/Mi are fixed, and then the 
value of mo is adjusted to give Odm^^ = 0.11. The four representative model lines were chosen 
to have (C24,-^o/Mi) = (0.19,1), (0.21,1), (0.21,1.5), and (0.24,1.5). The tuning of mo needed 
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FIG. 6: The branching ratios of scalar stoponium into the gg final state, for a variety of models of the 
type described in the text. In the left panel, —Aq/Mi = 1 at Mgut, with various C24 — 0.19 to 0.27 as 
labeled. In the right panel, C24 = 0.21, with various —Aq/Mi = 0.4 to 2.0 as labeled on the far right. In 
all cases, tan/3 = 10, ^ > 0, and Mi varies, with mo adjusted to give Quuh^ — 0.11. The models for which 
stop-mediated annihilations NiNi — > tt contribute more than 50% to l/flj^Mh^ are denoted by solid (blue) 
lines, other models by dashed (red) lines. The integrated luminosity needed to achieve a given Sj^pB can 
be obtained from fig. [3] by scaling by 1/[BR(77( 5<?)]^- 



is particularly mild in the regions indicated by the thicker solid (blue) lines for BR(r/^ gg), 
corresponding to models for which NiNi it dominates the annihilation of dark matter in the 
early universe. This typically gives m^, in the range 400-600 GeV. The model lines in fig. [5] show 
the common features that BR(?7^ gg) is quite high, typically 80% or higher for the range shown, 
and often in excess of 90% for smaller nirj- and smaller top-squark mixing. In the mass range where 
neutralino annihilation dominates, the branching ratios for W^W~ , ZZ, and h^h^ final states have 
a welcome dip, due to destructive interferences in the amplitudes for each of these final states. The 
dominance of the gg final state yields branching ratio to photons that are fairly constant, between 
0.004 and 0.005 over the relevant range of stoponium masses. Here and in the plots to follow, the 
lower endpoint on the model lines is set by the CERN LEP2 constraint on the Higgs mass, taken 
here to be m/jO > 113 GeV due to the theoretical errors in the computation. 

The crucial branching ratio into the gg final state is shown in fig. [6] for different slices through 
parameter space. In the left panel, —Aq/Mi = 1 is held fixed and C24 is varied over the range 
0.19 to 0.27 for which the stop-mediated neutralino annihilation to top quarks mechanism works 
efficiently, with mo always adjusted to give the observed dark matter density, and tan /? = 10 and 
/i > 0. Here BR(r/£ — > gg) is always high, but for larger C24, the stoponium mass is forced up by 
the LEP Higgs mass bound and will be difficult or impossible to observe at LHC for C24 > 0.26. 

In the right panel of fig. [6l we instead fix C24 = 0.21, and vary —Aq/Mi in the range from 0.4 
to 2.0. The BR(ry^ — > gg) tends to decrease slightly for larger —Aq/Mi (as the top-squark mixing 
increases), leading to enhanced branching ratios into W~^W~ , ZZ, and h^h^. However, the most 
dramatic effect is seen in the —Aq/Mi = 0.6 case, where the BR(r/j — > gg) becomes extremely small, 
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FIG. 7: Total integrated luminosity yielding an expected S/Vb = 2,3,4,5 for Mj^ in a bin with \Mj^ — 
m^jl < 0.02r7i^j, as a function of the stoponium mass m^^, for pp collisions at ^/s — 14 TeV, for compressed 
supersymmetry model lines of the type discussed in the text, with C24 — 0.21 and Aq/Mi = —1 (left panel) 
and An I My = -2 (right panel). 



due to the effect noted in ref. [9] of an s-channel resonance in stoponium annihilation to hh and ti 
from exchange. This can ruin the possibility of stoponium detection at the LHC. In the model 
framework with a common fixed mo and other values of C24, this resonant annihilation to quarks 
occurs for a range of —AqIM\ less than 1. It is interesting that these models correspond to the more 
optimistic projected sensitivity for the direct detection of dark matter in the next generation of 
low-background underground experiments 2j]. However, more generally, the heavy neutral Higgs 
boson mass can be made essentially arbitrary without changing the other essential features of the 
model, by assuming non- universal scalar masses at Mgut- Therefore, it is impossible to make any 
definitive statements about the complementarity of the detectability of the stoponium resonance 
at LHC and the direction detection of dark matter. 

Another dangerous decay mode is the /I'^/i'^ final state, which can dominate over all others if 
"m hPh^ is not too far above threshold, as noted in refs. [1, |l3]- In the dark-matter- motivated 
models I have studied here, this turns out never to be a fatal problem, because the stoponium mass 
is always sufficiently large. This is illustrated in fig. [TJ which shows the luminosity needed to obtain 
an expected S/\fB = 2, 3, 4, 5 for M^^ in a bin with |M^^ — m^^l < 0.02m^-, as a function of m^-, 
for two representative model lines with C24 = 0.21 and Aq/Mi = —1, —2. If the stoponium mass 
is small enough, detection might even occur with less than 10 fb~^ of data. In general, the models 
that are consistent with the dark matter scenario proposed in 23, l2j| can be prime candidates for 
stoponium detection in the diphoton mode, provided that the stoponium mass is not too large and 
not so close to m^jo as to allow a near-resonant annihilation decay. 
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FIG. 8: The branching ratios of scalar stoponium into gg, 77, W'^W~ , ZZ, and final states, for model 
lines motivated by electroweak-scale baryogenesis, as described in the text, with varying mf_^ . The left panel 
depicts a relatively optimistic case with m^o = 115 GeV, \Xt\/mf^ = 0.3 and the right panel a pessimistic 
case with m^o = 120 GeV, \Xt\/mf^ = 0.5. The range that can lead to electroweak-scale baryogenesis in 
the MSSM includes roughly 235 GeV < m^^_ < 270 GeV. 



IV. RESULTS FOR MODELS WITH ELECTROWEAK-SCALE BARYOGENESIS 

Another motivation for a relatively light top squark is the possibility of achieving electroweak- 



scale baryogenesis in the MSSM 2g]-[2g]. The necessity of a strongly first-order phase transition to 
a meta-stable electroweak symmetry-breaking vacuum limits the allowed parameter space, requiring 



a mostly right-handed top squark with mass less than nit. (For more details, see 26|- 28|.) Here, 



I will consider a model framework proposed in 27|], with an off-diagonal top-squark squared mass 



matrix element mtXt with 0.3 < \Xt\/mi^ < 0.5, and m^^ very large (here 10 TeV), tan/3 = 5 
to 10, m/jO between 115 and 120 GeV, and all other superpartners except the LSP supposed to 
be sufficiently heavy that they do not mediate large contributions to the stoponium decay width. 
Then the parameters with the most important impact on the stoponium decay widths are m^o, 
m^^, and the top-squark mixing angle. The region of parameter space where electroweak-scale 
baryogenesis can work is roughly 120 GeV < m^^ < 135 GeV 28|], but I will consider a wider range 



consistent with a meta-stable vacuum and the Higgs mass constraint from LEP2 ^27i |. 

In fig. [HI I show the relevant branching ratios for stoponium decay in a relatively optimistic 
case with m^o = 120 GeV and \Xt\/mi^ = 0.3 (left panel), and a pessimistic case with rrijyp = 115 
GeV and \Xt\/mi^ = 0.5 (right panel). In both cases, tan/3 = 10. Unfortunately, the branching 
ratio for the decay rji —>■ h^h^ is seen to be quite large above threshold [3], due to a small 
denominator (coming from the top-squark propagator) in the last term in eq. ()A.4p . The BR(r/^ — > 
hPhP) decreases as one moves to higher stoponium masses. I have made the optimistic but not 
unreasonable assumption that the other neutral Higgs boson is sufficiently heavy that the decay 
rjl — > hb is not near resonance and can be neglected. I have also optimistically assumed that the 
LSP is close enough in mass to ti so that r/^ NiNi is unimportant due to kinematic suppression. 
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FIG. 9: Total integrated luminosity yielding an expected S/vB = 2,3,4,5 for M^^ in a bin with \M^-y — 
m^jl < 0.02m,,j, as a function of the stoponium mass Wr/j, for pp collisions at ^/s — 14 TeV, for the two 
model lines depicted in fig. [S] The range that can lead to electroweak-scale baryogenesis in the MSSM 
includes roughly 235 GeV < m^, < 270 GeV. 



The existing collider lower limits on m^^ from the Tevatron and LEP2 do not constrain the top 
squark when m^^ —nij^^ is small, because of the softness of the charm jets from the decay. Obtaining 
a thermal dark matter density in agreement with WMAP from stop-neutralino co-annihilations in 
these models requires m^^ — m^^ ^20 GeV, but a smaller mass difference is still allowed, since it 
would just mean that the dark matter is something else, for example axions, or LSPs arising from 
a non-thermal source. 

Despite the large branching ratio to h?h? , there are still good prospects for stoponium detection 
at the LHC in these models, because the stoponium mass is necessarily not too large. This can be 
seen in fig. [9l which depicts the required luminosity for detection at various expected significances, 
i.e. the results of fig. H] divided by BR(77j gg)]^- In the most optimistic case, a clear observation 
of stoponium could be possible with less than 100 fb~^ over a large range of stoponium masses 
including the range that can accommodate electroweak-scale baryogenesis. If the decay h^h^ 
is kinematically forbidden, an observation could even be made with less than 10 fb^^. 

The detection of direct open light top squark production will probably be quite difficult at the 
LHC, especially if the m^^ ~ '^tvi ^^-ss difference is very small, because then the charm jets from 
ti cNi will be very soft. Most top squarks will likely come from gluino pair production followed 
by the decay g ^ tti. The Majorana nature of the gluino implies that half of the resulting events 
will have like-charge top quarks and soft charm jets. This was found to be a viable signal in \Q4\ . 
and endpoint analyses will allow the determination of relations between the gluino, lighter stop, 
and LSP masses. A relatively precise measurement of m^, would clearly be very helpful in pinning 
down the masses of all three particle. 
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V. OUTLOOK 

In this paper, I have examined prospects for observing stoponium at the LHC in the process 
^ ryj ^ 77, as first suggested by Drees and Nojiri in 0, 0]. As illustrated in two distinct 
motivated scenarios, this search will likely be a long-term project, requiring good precision in the 
electromagnetic calorimetry and high integrated luminosity. For a positive detection with 100 fb^^ 
or less, it is estimated to be probably necessary (but certainly not sufficient) that the stoponium 
mass is less than about 500 GeV, which happens to be in the middle of the range preferred by 



compressed supersymmetry models of the type discussed in refs. [23|, |2j]. If the stoponium mass 
is less than 300 GeV, detection might even be possible with 10 fb~^ or less. My estimates of 
the detectability for the idealized case of models for which BR(ry^ gg) w 1 are somewhat more 
pessimistic than in but there is a clear opportunity if Nature is kind. 

It would be useful to understand the higher-order corrections to stoponium production and 
decay, which might well have a large impact on the viability of the signal. Because these higher- 
order corrections are lacking for the signal, I have not included them for the background either, 
but they are likely to be substantial in both cases. In practice, the background will be obtainable 
from data using a sideband analysis. There is also an unknown, but possibly large, benefit from 
production of the excited states of stoponium adding to the signal. 

An important remaining question is how accurately the observation of stoponium can determine 
the lighter top-squark mass, and from it other superpartner masses. The answer is quite sensitive 
to unknowns, including the true size of the backgrounds, the experimental mass resolution and 
systematic errors for diphotons, as well as an estimate of theoretical errors in the binding energies 
and wavefunctions for stoponium bound states. However, it seems clear that the observation of 
stoponium would present a unique opportunity to gain precise information about the superpartner 
masses, and a key ingredient in deciphering the mechanism behind supersymmetry breaking. 



Appendix: Stoponium partial decay widths 

In this Appendix, I collect the results for the scalar stoponium annihilation decay widths. I 
agree with the results of ref. except for the gg and 77 widths; equations (A.l) and (A. 2) of 
that reference should each be multiplied by a factor of 1/2 on the right side. (This is on top of, 
and distinct from, the factor of 1/2 for identical particles which is correctly included in equation 



5) of ref. [9||.) For these two widths, I agree with refs. 11[ and 65|. Also, I have generalized ref. 

slightly, by including the effects of sbottom mixing and possible CP violating phases. I use a 
convention in which chargino and neutralino masses are always real and positive. All equations 
below use couplings and other notations and conventions as given in detail in the second sections 
of the two papers in ref. [66], which will not be repeated here for the sake of brevity. 
The general form for stoponium annihilation decay widths is 

for AB = gg, 77, ZZ, h^hP , Z7, W~^W~ , tt, bb, NiNj, with 6ab = 1 for the first four cases and the 
last case when i = j, and 5ab = in the others. Here A(x, y, z) = + ip + — 2xy — 2xz — 2yz, 
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and -R(O) is the wavefunction at the origin, for which I use the = 300 MeV parameterization 



from Table A.l of ref. 29[] in numerical work; see eq. (|2.3p of the present paper. It remains to give 
the spin-summed squared matrix element, 

For the gluon-gluon and photon-photon final states, the results are independent of soft 
supersymmetry-breaking parameters, by gauge invariance: 
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where qt = 2/2>. For the hPhP final state 
\M{r]i ^ h^h^)\^ = (A 



^<j,oiii*^hOhO<l,o 



j=l,2 •■■ii 



^-^ mi + mj — m 



(A.2) 
(A.3) 



(A.4) 



For the ZZ and W^W~ final states, 



ZZ) 



L \2 



2{a 



WW) 



(A.5) 
(A.6) 



where 
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i=i,2 



^2 + ^/2 



25 



/2 



4/2 2 
2 I 2 2 ' 

m^ -h m^ — m% 

ti tj z 



T 

"■WW 



ctww 



<^o=,o,^o 4m| - mjo ' 
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For the Zj final state, 

Y.\M{v,^Zj)f = 2qy{g'+g''){\L,f-Asl,/3r. 

The top-quark and bottom-quark final states have: 

\Mivi ^ tt)\^ = 6{ml - m?)(2ReK-] - mthtf + 24m|(ImM)2, 
^ \M{vi ^ bb)f = 6(m| - mg)(2ReK-,] - m,b,-,f + 24m| (Im[a,5])2, 



(A.7) 

(A.8) 
(A.9) 
(A.IO) 



(A.ll) 



(A.12) 
(A.13) 
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where 



att 



9 9 9 ~^ TT / 9 9 9 ; (A. 15) 

? + mi — mf 6 ^-^ mf + m — mf 





-I 2 ly - 12 I 



^_ ^ J- bOjt; ftOjfi f A 16) 



j=i ■ -c. "^h 

2 



3 ^2 _|_ j^2 _ j^2 

Finally, the final state with neutralinos has: 



(A.17) 



2m\{m\^ ^rn\^- Sm%m%^ - {m%^ + m%J/2 



+MibN^N^f]mfjmf^^{4.ml - 2vfi\^ - 2m\^ 
+2Re[o^^^^5^^^Jm^^(4m?^ - - 'im\^ 

+2Re[a^^^^5^^^Jm^^(4m| - SmJ^^ - vn\^, (A.18) 



where 



%Affc ~ ^2 _^™2_( 2 _^^2 W2 2^ 4^2 _7772 ' \ ■ ) 

N,Nk ^2 _^^2_('^2_ +^2_ I • J 

If there are no CP-violating phases, then CLfj^f^^ and bfj,fj^ are real, and eq. ()A.18P simplifies to: 
\M{vi ^ NjNk)? = 2[4m| - {m^^ + m^^f][aj^^^^ + (m^^. + m^j6^^.^^/2]2. (A.21) 
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